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ABSTRACT 



We here estimate the mass M of the central object in five Active Galactic Nuclei (AGNs) 
using the most recent reverberation data obtained by the AGN Watch consortium. The 
cross-correlation function (CCF) centroids of the broad Lya A1216 and C IV A1549 lines are 
used to estimate the size of the broad line region (BLR) in these sources. We calculate the 
velocity dispersions of these lines in the root mean square (rms) spectra, and then use our results 
to estimate M. We argue that our technique of calculating the velocity dispersion should work 
in the general case of an arbitrary line profile, unlike methods that depend on the measurement 
of the full-width at half-maximum (FWHM) of the broad line. We also show that our results 
agree with the FWHM method in the limit of a normal (Gaussian) line profile. The masses 
calculated here are considerably smaller than those calculated with the previous generation of 
reverberation data. 



Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — line: profiles 
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1. INTRODUCTION 

With our rapidly expanding knowledge of the Galactic Center (see, e.g., Melia 1999), we now have a 
real hope of unraveling how the BLRs in AGNs are produced, and how they evolve spatially and temporally. 
An understanding of the structure of the BLR in turn is important in determining how the gaseous 
environment behaves near the central engine, presumably a supermassive black hole (SMBH), in these 
sources. Unlike the Galactic Center, however, where we know the mass very accurately (e.g., Genzel, et al. 
1997; Ghez, et al. 1998), we do not yet have a reliable estimate of M for these objects. But since all of the 
relevant length scales depend on the central mass, one may reasonably expect that a careful modeling of 
these environments with temporally resolved line profiles can provide us with a tightly constrained range of 
values for M. 

Reverberation methods (e.g., Blandford & McKee 1982; Clavel et al. 1991) have become popular tools 
for estimating the BLR size in AGNs. The cross-correlation function (CCF) centroid T cent for a given line, 
typically measured in days, represents a delay between a change in the continuum flux and a corresponding 
(presumably causal) change in the line flux. From it, we can infer the responsivity- weighted mean radius, 
r ~ CT cent , of the line-emitting gas, assuming the delay is due to light-travel time effects (Koratkar & Gaskell 
1991a). If we make the very reasonable assumption that the central mass dominates the motion of the BLR 
gas, the velocity dispersion a 2 of this line can then be used to estimate the central mass, M w / cr ce „ t <7 2 / G , 
where / is a corrective factor of order unity representing our uncertainty in the geometry and isotropy of 
the BLR. 

An example of how this works in a simplified case, where the velocity distribution of the gas is such 
that the line profile in the root mean square (rms) spectrum assumes a roughly Gaussian shape, has 
recently been studied by Peterson & Wandel (1999). Using the UV and optical data for NGC 5548 from the 
extensive monitoring campaigns conducted by the AGN Watch consortium (see, e.g., Korista, et al. 1995; 
Peterson, et al. 1999), they computed the rms spectrum from the set of individual spectra. The line profiles 
in the rms spectrum represent those portions of the emission lines that are actually varying in response 
to the continuum fluctuations, and therefore correspond to the CCF measurements. Using this method, 
the FWHM of each line is measured directly from the line profile, and a velocity dispersion is calculated 
according to the relation a 2 = 3V^ WHM /4. Using the reported values of r cent , they obtained a central mass 
in this source of (6.8 ± 2.1) x 10 7 M Q . 

However, although estimating the velocity dispersion from the FWHM is valid in the case of a Gaussian 
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or other regular (e.g., Laplacian) line profile, this method cannot be used in the general case where the 
line profile is highly irregular. In this paper, we introduce a more generally applicable method to calculate 
the velocity dispersion directly, and show that our results agree with the FWHM procedure in the limit 
of a Gaussian profile. We avoid the simplifying assumptions concerning the line shape, which allows us to 
significantly broaden the class of objects for which the central mass may be estimated using reverberation 
methods. 

2. METHODOLOGY 

We have assembled and analyzed the UV data obtained with AGN Watch observations of five AGNs: 
3C 390.3 (O'Brien, et al. 1998), Fairall 9 (Rodriguez-Pascual, et al. 1997), NGC 3783 (Reichcrt, et al. 
1994), NGC 5548 during the HST FOS campaign (Korista, et al. 1995), and NGC 7469 (Wanders, et al. 
1997). The first object in this list is classified as a broad-line radio galaxy, whereas the others are Seyfert 1 
galaxies. The mean and rms spectra corresponding to each data set have been determined according to 
the procedure described in Rodriguez-Pascual, et al. (1997), and are shown for each of the five sources 
individually in Figures 1 and 2. These spectra have already appeared in the references cited above, but we 
reproduce them here as a class to demonstrate the range of line profiles one typically encounters with this 
process, and also to show the wavelength band used for each of the lines in our analysis (see below). We 
have selected the Lya A1216 and C IV A1549 lines for our analysis because they typically have the strongest 
signal-to-noise ratios. In those cases where the (redshiftcd) Lya line has been contaminated by geocoronal 
emission, we have attempted to set the lower bound of the wavelength bins beyond the contamination 
region. However, the Lya lines in NGC 5548 and NGC 7469 still show significant contamination (see 
Fig. 2). In all these objects, the Lya line is further contaminated by N V A1240, but no attempt to 
dcblend these two lines has been made here. The net effect of this is to increase our estimate of the Lya 
velocity dispersion. These contaminants make the Lya measurements much less reliable than the C IV 
measurements, which are only slightly affected by, e.g., He II A1640. 

These broad line profiles are understood to be produced by the composite emissivity of a broad line 
region (BLR) comprised of distinct clouds, which constitute a distribution in velocity around the central 
mass concentration. Along any given line-of-sight, say the z-axis, the velocity dispersion of these emitters is 
defined as 

a 2 z =< v\ > - < v z > 2 , (1) 
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where clearly v z denotes the z-component of the total velocity vector v. Then, assuming that \v\ «c, we 
can approximate the Doppler-shift formula as (A/A ) ~ l + v z /c, where Ao is the wavelength in the emitter's 
frame. Thus, as long as the velocity field is isotropic (i.e., a 2 = 3a 2 ), this relation can be used to replace v z 
in Equation (f ), and thereby to obtain an expression for a 2 in terms of observable quantities: 

a 2 = % « A 2 > - < A > 2 ) . (2) 

The mean values of A 2 and A appearing on the right hand side of this expression are often calculated from 
the data assuming that the line profiles are Gaussian, for which a determination of the FWHM is then 
sufficient. To keep this procedure as general as possible, we here instead calculate the expectation (or 
mean) value < g > of any given function g(X) of A using a flux- weighted sum: 

<g> = E a g(A»)Jline(A a ) ^ 
Ei^line(Ai) 

where Fi; no (A) is the (continuum-subtracted) line flux at wavelength A and the summation is taken over the 
appropriate wavelength bins Xi . The uncertainty in the calculated value of a is then estimated by assigning 
the highest and lowest plausible values for the underlying continuum. 

In our calculation of M , we have used the values of T cen t published in the references cited above for the 5 
sources considered here. The uncertainty in the value of T cen t is estimated using the Monte Carlo techniques 
described in White & Peterson (1994); however, rather than using a thin shell model for the BLR, we have 
instead assumed an extended BLR consistent with the observations and the general inferences drawn from 
modeling (see, e.g., Baldwin 1997). We chose the line response function to be ^(r, 9) = (1 + Acos6*)sin6l 
over the range of radii R in < r < R out , where A is the cloud emission anisotropy and Ri n and R out are the 
inner and outer radial extents, respectively, of the BLR. For example, a ^>(r, 9) that is fairly constant over 
all BLR radii is consistent with the density and velocity fields required to model the C IV A1549 line in 
Fairall 9 (Fromerth & Melia 1999). The inner extent of the BLR is set to correspond with the He II A1640 
cross-correlation lag (typically having the shortest lag of any of the broad emission lines), and the outer 
extent is optimized such that a histogram of modeled line delays peaks at the observed line lag (see White 
& Peterson 1994). We used the values of A ~ 1 for Lya and A — 0.7 for C IV A1549 calculated by Ferland, 
et al. (1992). 

Table 1 summarizes the results of these calculations. Our estimates of the errors associated with T cen i 
are much higher than those obtained using a thin-shell BLR model. This is intuitively correct, since we 
would expect an extended BLR to show a much broader, flatter CCF. In most cases, our extended BLR 
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model also yields larger error estimates for r cent than the model-independent method described in Peterson, 
ct al. (1998), which estimates errors based on uncertainties due to the observational sampling and the flux 
measurements only. This suggests that the latter method underestimates the actual uncertainty if the BLR 
is indeed spatially extended. The sources of error in a include the uncertain fitting of the continuum level, 
the blending of the broad line wings with the continuum, and line contamination (as discussed above). Since 
only the former is included in our error estimates, the o errors we quote here are probably underestimated. 
We also restate that the values of M listed in Table 1 are subject to corrective factors of order unity, 
stemming from our assumption of an isotropic velocity field. 

The last column in Table 1 lists the masses estimated by Wandel, et al. (1999) for the same objects 
using H/3 reverberation data. Their results assume a Gaussian line profile (usually a good approximation 
for the H/3 lines) and use the FWHM to estimate the velocity dispersion. They use the model-independent 
method of Peterson, ct al. (1998) to estimate the uncertainty in r ce „ t , so their reported mass errors may be 
underestimated (as discussed above). 

3. ANALYSIS & CONCLUSIONS 

Perhaps the most striking inference we can draw from an inspection of Table 1 is that the masses 
calculated here are considerably lower than those estimated using older reverberation data. For example, 
Koratkar & Gaskell (1991b) predicted masses of x 1O 7 M , x 10 s M Q , and 4.7t^ x 10 8 M Q 

for NGC 3783, NGC 5548, and Fairall 9, respectively. The main discrepancy seems to be that the higher 
(and regularly-spaced) sampling rate adopted by the AGN Watch has facilitated the detection of shorter 
time scale variations and cross-correlation lags for these objects. An additional effect stems from our use 
of the rms spectrum rather than the mean spectrum to estimate T^fwhm or a. A cursory glance at Figures 
1 and 2 suggests that the velocity dispersions of the rms spectral lines are smaller than those of the mean 
spectra. The physical interpretation of this effect is not yet clear, though it may be signaling that the 
emission from the innermost BLR gas does not vary strongly in response to continuum changes. 

Our mass estimates agree with the Wandel, et al. (1999) results within error bars in all cases except 
the Lya results of 3C 390.3 and NGC 3783. The NGC 5548 mass estimate is also consistent with the result 
obtained by Peterson & Wandel (1999) using the same spectral data, though with the assumption that the 
lines are Gaussian. Collier et al. (1998) estimated the mass of NGC 7469 using the spectral variability in the 
optical regime combined with a measurement of Vfwhm for the Ha and H/3 lines to obtain a mass for this 
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object of 8.3 ± 2.4 x 10 6 M©, which is again consistent with our result. These cases provide some support for 
the validity of our method, since they show a consistent mass determination for sources that have "normal" 
(i.e., Gaussian) line profiles. Our method is clearly more general, since it appears to be applicable also to 
sources in which the velocity distribution of the BLR does not naturally produce a Gaussian line shape. 

A caveat one must consider with this procedure is that the error associated with the Lya fit can be 
significant. The mass estimates made using the Lya data are higher than those based on the C IV line in 
all cases except Fairall 9. As we discussed above, this is probably due to a contamination of the line, which 
suggests that the C IV results may be more reliable. It is possible that this problem may be alleviated by 
line deblending techniques, which we have not incorporated into our analysis here. 

Our results support the argument made by Peterson & Wandel (1999) that a very large mass must be 
present within a distance of a few light-days from the centers of many AGNs. For the class of 5 objects 
in our sample, this mass is at least 8 x 10 6 M©, which is rather compelling evidence for the presence 
of a point-like object. This mass density is even higher than that in the Galactic Center, where we are 
reasonably sure that the enclosed mass within 0.015 pc of the nucleus is close to 2.6 x 10 6 M©. Alternative 
models for the latter that invoke a distribution of dark objects, such as neutron stars, white dwarfs, or 
~ 10 M© black holes, can now be ruled out on the basis that in equilibrium, such a distribution would 
have a core radius significantly larger than the observed limiting size of 0.015 pc. In the sample of AGNs 
considered here, this exclusion of a distributed dark matter contribution to the enclosed mass appears to 
be even stronger (cf. Maoz 1998). 

One of the most important results presented here is the relatively accurate determination of the central 
point mass in these objects. With the exclusion of 3C 390.3, which appears to belong to a different mass 
class, the masses of the other four AGNs (all Seyfert I's) are very close to each other and are near the low 
end of the mass scale for AGNs. This characteristic, and the fact that the SMBH at the Galactic Center has 
a mass within a factor of 4 or so of this range, hints at a possibly similar formation and evolutionary history 
for these nuclei. In the Galactic Center, the power generated by the supermassive black hole seems to be the 
quasi-steady accretion of plasma ejected by wind-producing stars in the circumnuclear environment, with 
a rate M gc « 2 x 10~ 4 M© yr" 1 (Coker & Melia 1997). Thus, given a central mass M gc w 2.6 x 10 6 M©, 
the characteristic time M gc /M gc w 13 billion years suggests that the black hole Sgr A* may have grown 
to its current size gradually over a Hubble time via the steady accretion of gas from its environment. We 
speculate that the nuclei in the 4 Seyfert I's we have considered here may themselves have formed and 
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evolved in a similar fashion, and that the range in masses may therefore be due principally to a variation of 
the nuclear gas density from which the black holes accrete. 
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Table 1. Velocity Dispersions and Mass Estimates of AGNs. 







Tcent 




Mass b 




M w c 


AGN 


Emission Line 


(days) 


(km s- 1 ) 


(10 7 M Q ) 


Reference 


(1O 7 M ) 


3C 390.3 


Lya A1240 + N v A1240 




TMOtgg 


70±H 


O'Brien et al. 1998 


391^1 




Civ A1549 


4ll* 4 


8iio±S8 


53l? 3 






Fairall 9 


Lya A1240 + N v A1240 




6180±||8 


11+8 


Rodriguez-Pascual et al. 1997 


8.7ti$ 




Civ A1549 


32±«t 


43801™, 


1 2 +io 

1Z -10 






NGC 3783 


Lya A1240 + N v A1240 


4.71? 


4140l 6 9 3 ° 


1 6 +0 ' 8 

i -°-0.4 


Reichert et al. 1994 


i-Hi-.o 




Civ A1549 


3.9±\ 




-i O+0.8 






NGC 5548 


Lya A1240 + N v A1240 


6.9l 3 2 


65201J™* 


5 7+3-1 

°-'-2.4 


Korista et al. 1995 


6-s±i:g 




Civ A1549 


7.0l 3 2 


63201^8 


c c+3.3 
o.o_2.g 






NGC 7469 


Lya A1240 + N v A1240 


2.31} 


54601$* 


i-3±8:S 


Wanders et al. 1997 


0.7618;^ 




Civ A1549 


2.71} 


3970l«8 


o.83i8:S 







a The errors in a are probably underestimated due to wavelength bin assignment and line contamination. 

b Mass estimates are subject to systematic error due to our assumption of isotropy. 

c The mass estimates reported in Wandel, et al. 1999 based on H/3 reverberation data. 

tRodriguez-Pascual, et al. 1997 attach no significance to this value. 

*The Lya profile in this source is badly contaminated by geocoronal emission. 



- 12 - 



i — n—p — 1 — i — 1 — 1 — 1 — r-fr 1 — 1 — 1 — r 

A| 3C 390.3 




j i i i i i i i i i i i i 



1200 1400 1600 1800 

Wavelength (A) 



Fig. 1.— Mean (light plot) and rms spectra (dark plot) of 3C 390.3, Fairall 9, and NGC 3783. The line 
flux wavelength ranges for Lya + Ny and Civ arc indicated. Note that, in general, the lines are not well 
described by Gaussian. 
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Fig. 2.— Mean (light plot) and rms spectra (dark plot) of NGC 5548 and NGC 7469. The line flux 
wavelength ranges for Lya + Ny and Civ ar e indicated. Note that the Lya line is show significant galactic 
contamination in both cases. 



